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Summary. The expansion of the electron-positron plasma in the GRB phenomenon 
is compared and contrasted in the treatments of Meszaros, Laguna and Rees, of 
Shemi, Piran and Narayan, and of Ruffini et al. The role of the correct numerical 
integration of the hydrodynamical equations, as well as of the rate equation for the 
electron-positron plasma loaded with a baryonic mass, are outlined and confronted 
for crucial differences. 



1 Introduction 

One of the earliest contributions in the theoretical understanding of GRBs is 
often referred to the work of [1]. Some crucial ideas which later became very 
important in this field were there presented (e.g. the relevance of an electron- 
positron plasma, the possibility of baryonic loading of such plasma, and the 
evolution from an optically thick to an optically thin phase). The treatment, 
based on qualitative consideration, provided little information about the dy- 
namics. The electron-positron plasma was there purported to evolve toward 
a proliferation of photons and pairs, leading to a degradation of the mean en- 
ergy of particles in the plasma and a consequent cooling. Such a plasma would 
become transparent on a very short time scale and no dynamical phases are 
envisaged making this model inappropriate for GRBs. The physical reasons 
of the non-validity of these crucial conclusions are being discussed elsewhere 
[2]- 

It soon became clear, however, following the work of [3] and [4] that the 
major characteristic of a sudden energy release process in electron-positron 
plasma leads to a very rapid self acceleration of a shell of material, reaching 
ultra-relativistic regimes with Lorentz gamma factors in the range 10 2 < 7 < 
10 3 . The major results were obtained in [5], in [6], in [7] and in [8, 9]. 

In our model the analysis of the dynamical expansion of the electron- 
positron plasma is not just a topic of academic interest, it is indeed crucial 
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to the description of the entire GRB phenomena. The initial process is the 
vacuum-polarization around the black hole [10], followed by the dynamical 
expansion of the pair plasma, leading to an ultrarelativistic accelerated motion 
[8, 9]. Such acceleration leads to a shell of baryons with Lorentz gamma factor 
7 ~ 100 — 300 which, by interacting with the interstellar matter (ISM), gives 
origin to the afterglow. A first crucial signature in our model is carried by 
the radiation emitted at the moment of transparency, what we have called 
the Proper-GRB. The current models in the literature identify two different 
components in the long GRBs, the "prompt radiation" and the "afterglow", 
the first one being emitted by an "inner engine" (see e.g. [11] and references 
therein). In our model, instead, the "prompt radiation" is an essential part 
of the "afterglow" and the entire long GRB phenomena is uniquely due to 
external shock processes. 

For the above reasons we return, in this communication, to a comparison 
and contrast between our results and the ones in the current literature. Now 
we can address with clarity the analogies and the differences in the treatments 
and pin down the source of such differences as due to inaccurate theoretical 
work or / and equally inaccurate numerical computations. 

2 Contrasts in the dynamical description of the 
expanding plasma 
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Fig. 1. Left: Relative energy release in the form of photons emitted at transparency 
point _E 7 (solid line) and kinetic energy of the plasma Ek (dashed line) of the baryons 
in terms of the initial energy Eq of the electron-positron plasma as computed in our 
quasi-analytic model [8, 9]. The quantities are given as a function of the baryon 
loading parameter B. The bold line denotes the total energy of the system in terms 
of initial energy Eq which is, as it should, constant and equal to 1. Right: The same 
quantities are computed in the Meszaros, Laguna and Rees model [5]. Note that in 
Ref. [5] the parameterization is done as a function of the dimensionless parameter 
rj = 1/B. The bold line denotes the total energy of the system in terms of initial 
energy Eq, which is not conserved in such a treatment [5]. 
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Fig. 2. Left: Rclativistic gamma factor when transparency is reached as a function 
of the baryon loading parameter B. The bold line denotes the numerical results 
obtained in Ref. [8, 9], taking into due account the electron-positron rate equation. 
The plain line corresponds to the analytical estimate from Shemi and Piran model 
[12], neglecting the rate equation. The dashed line denotes the asymptotic value of 
the Lorentz gamma factor at transparency 'jasym = B~ (see Ref. [9] for details). 
Right: Relative energy release in the form of photons emitted at transparency point 
_E 7 of the GRB in terms of initial total energy Eq depending on the baryon loading 
parameter B. The bold line represents our numerical results, already given in Fig. 
1-Left. The plain line shows the results for the analytic model of Shemi and Piran 
[12]. The dashed line shows the difference between our numerical analysis, taking 
into proper account the rate equation [8, 9] and the approximate analytical [12] 
results. 



Our results are represented in Fig. 1-Left, showing the energy release at 
transparency carried by photons and in the form of kinetic energy of the pulse, 
depending on the loading parameter B = Mc 2 /Eq, where M is the total mass 
of the plasma, c is the speed of light, Eq is initial energy of the system. In 
Fig. 1-Right we represent the corresponding results by [5] which are markedly 
different from ours, and do not fulfill the basic requirement of the conserva- 
tion of energy. In Fig. 2-Left we represent the results of [12] about the gamma 
factor at transparency which overlap with ours. The major difference between 
our treatment and the one by [6] and collaborators stems from different ac- 
curacy in the description of the electron-positron pairs equations which we 
have explicitly integrated in full details with the dynamical equations of the 
pulse. In Fig. 2-Right we represent some analogies and differences between 
the two treatments. In particular, the energy release in the form of photons at 
transparency, which is crucial for the detection of the P-GRB, is qualitatively 
correct in Ref. [12] but underestimated in view of the simplified analytical 
approach adopted, which does not take into proper account the explicit inte- 
gration of the electron-positron rate equation. The moment when the pulse 
reaches transparency and the corresponding radius are also different in simple 
analytical models and our detailed numerical computations, which leads to 
different predictions for the energetics. 

In addition to the above differences, the treatment in [9] clearly predicts 
an instability in the expanding plasma for a value B > 10~ 2 . 
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3 Conclusions 

The analysis of the electron positron plasma expansion, far from being purely 
academical, is essential in describing the entire GRB phenomenon in our 
model. We are currently working in proving the uniqueness of our model. 
It is clear that any minimal deviation in the integration of the hydrodynami- 
cal equations, or any inadequacy on the integration of the rate equation, may 
lead to the impossibility of reaching the correct theoretical model of GRBs by 
lacking the necessary accuracy in the description of the fundamental process 
determining the energetics of GRBs. Indeed the agreement of our theoretical 
model [13, 14, 15, 16, 17, 18, 19] with the observations has been successfully 
tested in four different sources for the intensities of both the P-GRB and the 
afterglow in selected energy bands. These sources are: GRB 991216 [20], GRB 
980425 [21], GRB 030329 [22], GRB 031203 [23] (see Fig. 3). It is particularly 
interesting that, in all the above sources, the value of the B parameter is 
smaller than 10~ 2 as clearly predicted in Ref. [9]. 
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Fig. 3. The relative intensity of the P-GRB and the afterglow are given for four 
selected sources as a function of the baryon loading parameter B. 
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